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A description of fragmentation functions which satisfy the momentum and isospin sum rules is 
presented in an effective quark theory. Concentrating on the pion fragmentation function, we first 
explain why the elementary (lowest order) fragmentation process q — > qn is completely inadequate to 
describe the empirical data, although the "crossed" process -n — » qq describes the quark distribution 
functions in the pion reasonably well. Taking into account cascade-like processes in a generalized jet- 
model approach, we then show that the momentum and isospin sum rules can be satisfied naturally, 
without the introduction of ad hoc parameters. We present results for the Nambu-Jona-Lasinio 
(NJL) model in the invariant mass regularization scheme and compare them with the empirical 
parametrizations. We argue that the NJL-jet model, developed herein, provides a useful framework 
with which to calculate the fragmentation functions in an effective chiral quark theory. 
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I. INTRODUCTION 



Quark distribution and fragmentation functions are the basic nonperturbative ingredients for a QCD-based anal- 
ysis of hard scattering processes [H, 0, H, 13, IE @|- Distribution functions can be extracted by analyzing inclusive 
processes @, H| and their description in terms of effective quark theories of QCD has been quite successful [H, [l(| • In 
recent years there has been a significant effort to extract the fragmentation functions by analyzing inclusive hadron 
production (semi-inclusive) processes in e + e~ annihilation, deep-inelastic lepton-nucleon scattering and proton-proton 
collisions [111 H^ |. Besides being of fundamental interest in their own right, knowledge of fragmentation functions is 
essential for the extraction of the transversity quark distribution functions @, [l3| from data, and to analyse several 
other interesting effects in semi-inclusive processes [14| . 

Because of the importance of the fragmentation functions many attempts have been made to describe them using 
effective quark theories 15]. However, in order to achieve reasonable agreement with the empirical parametrizations 
it was necessary to introduce new parameters, like normalization constants, which cannot be justified on theoretical 
grounds. A description of fragmentation functions within effective quark theories, which automatically satisfies the 
relevant sum rules [H and describes the empirical data reasonably well - without introducing new parameters into 
the theory - has hitherto not been achieved. 

This failure to describe the fragmentation functions in the same framework which is successful at describing the 
distribution functions is surprising, because there exists a general relation, the so called Drell-Levy-Yan (DLY) relation 
[lfjl [l7j , which suggests a way to compute the fragmentation functions by analytic continuation of the distribution 
functions into the region of Bjorken x > 1. Although the derivation of this relation appears to be very general (as 
we show in Appendix A), the basic assumption that the distributions and fragmentations are essentially one and the 
same function, defined in different regions of the scaling variable, has not been proven. Moreover, the approximations 
used to calculate the distribution functions may not be sensible for the fragmentation functions and vice versa. For 
example, in the fragmentation process of a quark into a pion, q — > n + n, where n is a spectator, there is no a 
priori reason to truncate n to a single quark state, as the DLY crossing arguments would suggest for the case of 
a Bethe-Salpeter type vertex function for n — > qq. One can actually give a quantitative argument that the lightest 
component of n is dominant only if the scaling variable z is very close to unity flU . 

On the other hand, the phenomenological quark jet-model, as formulated originally by Field and Feynman [l9l |. 
suggests that the meson observed in a semi-inclusive process is one among many, that is, the spectator state n contains 
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many mesons. This model is based on a product ansatz for a chain of elementary fragmentation processes, where 
in each step a certain fraction of the quark momentum is transferred to a meson, until eventually a very soft quark 
remains. This final quark is assumed to annihilate with the other remnants of the process without producing further 
observable mesons. 1 In order for all of the quark light-cone momentum to be transferred to the mesons, it is actually 
necessary to assume an infinite number of steps (mesons) in the decay chain, as will be explained in more detail 
in Section llVI In this case, it is possible to satisfy the momentum sum rule for fragmentation functions [3|, which 
is assumed valid in QCD-based fits to the data [ill, EH- Clearly, this sum rule cannot be satisfied in a single step 
elementary fragmentation process. 

The purpose of this paper is to apply the method of the quark jet-model to calculate the spin-independent frag- 
mentation functions in an effective chiral quark theory, which has proven to be very successful for the description of 
quark distribution functions I22I [23j . We will concentrate on quark fragmentation into pions within the Nambu- 
Jona-Lasinio (NJL) model [24j, however the methods illustrated here can easily be extended to other fragmentation 
channels and applied within other effective quark theories. In order to reconcile the quark jet-model with our present 
NJL model description, we will introduce a generalized product ansatz, which allows for the fragmentation of a quark 
into a finite number of pions according to a certain distribution function, and in the end we take the limit of infinitely 
many pions. We will show how the momentum and isospin sum rules emerge naturally without introducing any new 
parameters into the theory. Our numerical results will demonstrate that this NJL-jet model provides a very reasonable 
framework for describing the fragmentation functions. 

This paper is organized as follows: In Section [TT] we begin with the operator definitions for the quark distribution 
and fragmentation functions and move on to discuss the sum rules and the DLY relation. In Section IIIII we give the 
expressions for the elementary fragmentation functions in the NJL model and discuss their physical interpretations and 
sum rules. In Section HVl we introduce the generalized product ansatz to describe a chain of elementary fragmentation 
processes in the spirit of the quark jet-model, derive the integral equation for the fragmentation function and discuss 
the momentum and isospin sum rules. In Section [V] we explain the model framework for the numerical calculations, 
present results and compare them with the empirical fragmentation functions. A summary is given in Section [VII 

II. OPERATOR DEFINITIONS AND SUM RULES 

Operator definitions and sum rules for fragmentation functions were first given in Ref. !3J and were further elucidated 
in Ref. (25j . In this Section we summarize the basic relations for the fragmentation functions and for clarity include 
those for the distribution functions also. The spin-independent distribution function of a quark of flavour q inside a 
hadron of spin-flavour h (for example h = p|,7r + , etc.) and the spin-independent fragmentation function for q — * h 
are defined by 

(*) = \j ^^^^ ii ( P (^)l?(o)K) 7 + (pnMu-Mh)), (1) 

D " (Z) = T2 I ^ e4p - W " /Z E„(PW,PnR0)|0) 7 + (0|^^-)b(/ l ),^). (2) 

The field operators refer to a quark of flavour q, although it is not indicated explicitly. The symbol p(h) refers to a 
hadron h with momentum p and p n labels the spectator state. The light-cone components of a 4- vector are defined as 
= (a + , a~ ,ax) with = (a°±a 3 )/v / 2- Covariant normalization is used throughout this paper and the summation 
symbol J\ includes an integration over the on-shell momenta p n . 2 Both expressions in Eqs. ([I]) and ([2]) refer to a 
frame where p T — 0. The physical content of the functions in Eqs. |T]) and ^ is most transparent if we introduce 
the "good" light-cone quark field ip+ [22L I26L l27l ]. which is defined by ip+ = A+tp where A+ = ^7~7 + and can be 
expressed as the Fourier decomposition 

<M<0 = -7= / (|^7a E (M fc ) u +a {k) e- ik -»- + 4(fc) v +a (k) e^) . (3) 



1 This picture of independent fragmentation is appealing because of its simplicity. More elaborate models for hadronization are the string 
model [20| or the cluster model |2l| , which are suitable for Monte Carlo analysis. 

2 In this normalization (p' (h')\p(h)) = 2p_ (2-7r) 3 <5( 3 ' (p' — p) 5 hh i and \p(h),p n ) = \/^(^ n ) 3 P- a I(p)|Pi)> wrtn [ a h(p'), a h(p)]± = 
5( 3 ) (p' — p). The summation defined by = / pjjp -M^) ©(Pno), where M n is the invariant mass of n, can also 
be expressed in terms of light-cone variables. 
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The index a denotes the spin-color of a quark with flavour q and the spinors are normalized as u^ a i (k) u+ a (k) — 
v +u'(k) v +a{k) — y/2k- S a '. a- Substituting these expressions into Eqs. (Q} and ([2]) and using the result -07 + -0 = 
\Z2iJ)\_ip+, gives the following relations which are independent of the normalization of states: 

Here dx = dk-/p-, that is, k- = xp- for some fixed p_ > and dz = dp-/k—, implying p_ = zfc_ for some fixed 
fc_ > 0. The creation and annihilation operators, a h and at, refer to the hadron h (see footnote and k(a) labels a 
quark state of flavour q with momentum k and spin-color a. 

According to Eq. Q we can interpret fq(x) as the light-cone momentum distribution of g in h, where a sum 
over the spin-color of q is understood, while the spin of h is fixed. However, the result is independent of this spin 
direction, since we will only consider the spin-independent distributions. As mentioned earlier, Eq. ([5]) refers to the 
frame where the produced hadron h has p T = 0, but the fragmenting quark has non-zero fcy. To interpret this result 
as a distribution of h in q, it is necessary to make a Lorentz transformation to the frame where k± = 0, but h has 
non-zero p ± (note the distinction between the subscripts T and _L). This is discussed in detail in Refs. [H, with 
the result that one can simply substitute 

k T = ~ E ±, (6) 

z 

leaving everything else unchanged. We then obtain from Eq. ([5]) the result 

D h (z)dz = -dp / d 2 ( 7 ) 
u q[ z } az Q ap-ja P± 2_^ (k(a)\k(a)) ' () 

where the fragmenting quark now has k± = 0. According to Eq. (J7J| we can interpret D^(z) as the light-cone 
momentum distribution of h in q, where the factor 1/6 indicates an average [25j of the spin-color of q, while the spin 
of h is fixed. 3 In fact, for the elementary distribution and fragmentation functions considered in the next section, the 
naively expected relation 

D%(z) = j-f%(z), (8) 

is valid. Where dh is the spin degeneracy, or, in the general case, the spin-color degeneracy of h. Generally however, 
this relation is not necessarily valid, because q is off-shell (its virtuality being determined kinematically by the scaling 
variable and the transverse momentum) and h is on-shell, which breaks the naive symmetry under the interchange 
q h. 

To obtain the momentum sum rule from Eq. (0) we multiply both sides by z = p_/fc_, integrate over z from to 1 
and sum over h. 4 Then one notes that the momentum operator, represented in terms of hadron operators, is given by 

p- = J2 1 d p- 1 d V (p- 40) a h(p)) ■ (9) 



By assuming that the quark state \k(a)) in Eq. |J7J| is an eigenstate of this operator with eigenvalue we obtain the 
momentum sum rule 



^ [ dzzD'^(z) = 1. (10) 
Jo 



h 



3 For the generalized case where h can also be a quark, we summarize the definitions as follows: fq(x) refers to fixed flavours of q and 
h, while all other quantum numbers of q (spin, color, etc) are summed over, with those of h are fixed. D™(z) refers to fixed flavours 
of q and h, with an average over the other quantum numbers of q (spin, color, etc), while those of h are fixed. This definition has the 
advantage that in a semi-inclusive process, which involves the product f^(x)Dg(z), the quark spin-color summation is included in / 
but not in D, which avoids double counting. 

4 A subtle point here is that in order to get an integral dp- on the right hand side of Eq. J9}, one has to choose k— = oo. This does 
not influence the result, which depends only on z. 
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The physical content of Eq. (p~0|) is that 100% of the initial quark light-cone momentum is transferred to the 
hadrons. The condition which lies at the basis of Eq. (jTUJ) is that the initial quark state is an eigenstate of the 
momentum operator, Eq. ([9]), expressed solely in terms of hadrons. That is, the quark hadronizes completely in the 
sense that it gives all of its light-cone momentum to the hadrons. 
A similar argument leads to the isospin sum rule [3|, namely 

V f dzt h D h q {z) = t q: (11) 
h Jo 

where t q and th denote the 3-components of the isospins of q and h. The physical content of this sum rule is that all 
of the isospin of the initial quark is transferred to hadrons, which is possible since the definition in Eq. (|2|) implies an 
average over the isospin of the soft quark remainder of a fragmentation chain (see Section HV)) . In general, there is no 
sum rule for the baryon number or electric charge, because the baryon number or average electric charge of the quark 
remainder is not zero. If we simply integrate both sides of Eq. ([7]) over z, we get the hadron multiplicity, which can 
be interpreted as the number of mesons per quark. However, there is no conservation law which leads to a sum rule 
for the multiplicity. 

There is an interesting relation based on charge conjugation and crossing symmetry, between the fragmentation 
function for physical z < 1 and the distribution function for unphysical x > 1: 




which is called the Drell-Levy-Yan (DLY) relation [n| [ItJ • Here s q and Sh are the spins of q and h respectively, and 
d q is the spin-color degeneracy of q. We derive this relation using two independent methods in Appendix A. The first 
approach, which follows the original arguments [l6j| . compares the hadronic tensors for eh — ► e'X (inclusive DIS) and 
e + e~ — > hX (inclusive hadron production), and uses crossing relations for matrix elements of the current operator. 
The second method - which to the best of our knowledge has not been published before - starts directly from the 
operator definitions in Eqs. (fTJ) and ^ and uses charge conjugation and crossing symmetries for matrix elements of 
the quark field operator. If one has an effective quark theory to calculate the quark distribution functions, Eq. (|12|) 
would suggest a straightforward way to obtain the fragmentation functions. However, as will become clear in the 
following sections, for the lowest order (elementary) processes such an attempt leads to disastrous results. That is, the 
fragmentation functions obtained in this way are one or two orders of magnitude smaller than the empirical functions 
and the sum rules in Eqs. (p~0|) and (fTTj) are not satisfied. 

The reasons why Eq. (fT2|) fails in actual applications are as follows: (i) It is based on the assumption that the 
distribution functions can be continued analytically beyond x = 1. However, it is well known that the Q 2 evolution 
equations lead to singularities at x = 1, which are (regularized) infrared singularities arising from the vanishing gluon 
mass 0, Hll . These render an analytic continuation impossible. Someone may still argue that Eq. (jX2j) should be used 
only at the low energy (model) scale, however it is actually broken there also, because of the cut-off regularization. We 
will discuss this point in detail in the next section, (ii) Most importantly, approximations which work reasonably well 
for the distribution functions may not be sensible for the fragmentation functions and vice versa. For example, the 
assumption that the pion is a qq Bethe-Salpeter bound state is very reasonable for the distribution function 0, ficj ] , 
but the DLY crossing arguments then imply the truncation of the spectator state p n , of Eq. @, to a single quark 
state. Although this simple assumption does not lead to any violation of conservation laws, the sum rules in Eqs. (jTUJ) 
and (jlip cannot be satisfied in a single step fragmentation process. 

For these reasons, we will not rely on Eq. (fT2| to calculate the fragmentation functions, although we will confirm its 
formal validity for the lowest order (elementary) functions. We note that the arguments given above do not question 
the usefulness of Eq. p^|) as a means to relate the kernels of the Q 2 evolution equations for the distribution and 
fragmentation functions (see Ref. [l7| and Appendix B). In fact, it is known that at leading order (LO) in a s this 
relation between the kernels is valid, although it is violated at next-to-leading order (NLO) [291 ]. 



5 The average of the electric charge is zero only if SU(3) flavour symmetry is assumed. 
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FIG. 1: Figure (a) depicts the cut diagram (left) and Feynman diagram (right) for the distribution function fq(x). Solid lines 
denote the quark and dashed lines the pion. Here fc_ = ip_ and the two quark lines with momentum k are connected by a 
7 + . Figure (b) depicts the cut diagram for the fragmentation function dj(z). Here p_ = zk- and the two quark lines with 
momentum k are connected by a 7 + . This diagram refers to a frame where p T — and the substitution given in Eq. ((6]l is 
performed in the final transverse momentum integral. 



III. ELEMENTARY DISTRIBUTION AND FRAGMENTATION FUNCTIONS 



The elementary distribution and fragmentation functions for the pion are represented in Figs. [T] as cut diag rams. 
Since the distribution function can also be obtained from a straightforward Feynman diagram calculation |22l |30| , 6 
we also illustrate the Feynman diagram for the distribution function on the right hand side in Fig. [T^. We denote 
the elementary fragmentation function by d^ in order to distinguish it from the total fragmentation function 
determined in Section [TV] We obtain the following expressions from the diagrams in Figs.[T] 7 



d A k 



Ttd [5i.(%+5F(fe)75 {t-f- M) 75] S(k- - p-x) 5(( P - kf - M 2 ) 



~ (1 + r„r q ) 6gl J 



1 2 Z 



M 2 



( 27r ) 3 [fey + M 2 — m%x(l — a;)] ! 
d A k 



1 

x = — 

z 



6 Jq 

1 . , 2 

: - (1 + T 7r r q ) zg„ 



(27T) 



d 2 p± 



Tv D [S F (k) 7 +S F {k) l5 (k'-p:-M) l5 ] 8(k- -p-j z) S((p — k) 2 - M 2 ) 



Pi 



(2tt) 3 [pi + M 2 z 2 + {l-z) ml] 



(13) 
(14) 

(15) 
(16) 
(17) 



where Ttd indicates a trace over Dirac indices only. The Feynman propagator of a constituent quark with mass M 
is denoted by Sf and g v is the pion-quark coupling constant. In the NJL model g v is defined via the residue of the 
qq i-matrix at the pion pole, and can be expressed in terms of the qq bubble graph by 



dll^q 2 ) 



dq 2 



where II 7r (q 2 ) = 6i 



d A k 



Tr D [j5SF(k)j 5 S F (k + q)] . 



(18) 



We use the isospin notations (t U5 t^) = (1,-1) and (t w + , t v o , T n - ) = (1,0,-1). For the distribution function in the 
physical region (0 < x < 1) a factor 6 (p_ — fc_) = 9(1 — x) has to be supplied in Eq. (flljj) . which expresses the 



This is seen simply by using completeness in Eq. {T} and the identity i])+(fi)^i])+(ui~ ) = T (i/i-|-(0)t^>+(u — )) in the limit w+ — + — e, 
which follows from causality. 

The expressions given in this section refer to the NJL model, however they actually have the same form in any effective chiral quark 
model with point-like pion-quark vertex functions. 
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— > 

FIG. 2: The quark self-energy, Eg^fc) = — 3ig% J T^rz 'fsSF^k — p)7sAf(p), where Ap is the Feynman propagator of the pion. 



fact that the intermediate antiquark in Fig. [T^, has positive energy. Similarly, for the fragmentation function a factor 
0(/c_ —p~) = 0(1 — z) has to be supplied in Eq. (|15p . because the intermediate quark in Fig.[]J> has positive energy. 
To obtain Eq. (fT7j) we made the substitution given in Eq. ([6]) . 

The DLY relation on this level, indicated in brackets as Eq. ([T5|. shows that Eq. (|T5|) can be considered as a 
generalized distribution function, which gives the physical distribution function in the region < x < 1 and the 
fragmentation function in the region x = 1/z > 1. The reason why we indicate this relation only in brackets is that 
it is violated if the integrals are regularized. For example, if we use a sharp cut-off (A) for the transverse quark 
momentum in Eq. (114[) . a strict application of the DLY relation would mean that the transverse momentum of the 
produced pion in Eq. (|17|) should be cut at z A, which is unacceptable. The more physical procedure is to impose 
|fcr| < A on Eq. (fT4")l and < A on Eq. (fTT)) , which breaks the DLY relation. A similar breakdown of the DLY 
relation occurs in any other sensible regularization scheme. A noticeable consequence of this is that in the chiral limit 
the distribution function of Eq. (fT4]) becomes a constant, but the fragmentation function of Eq. (fTT)) is not linear in 
z, as the DLY relation indicated in Eq. (fT5|) would suggest. 

The relations for the distribution function 

f 1 1 f 1 1 1 

J dxfg(x) = - (1 + r^Tq) , and J dxxfg(x) = - (1 + T 7T T q ) • -, (19) 

lead to the usual number and momentum sum rules. For the elementary fragmentation function the following relation 
is obtained from Eq. (JTjJ : 

1 dzd%(z) = \(l + T„T v )(l-Z Q ) / dzJ2<r q (z) = l-Z Q , (20) 

where Zq is the residue of the quark propagator in the presence of the pion cloud. It is expressed in terms of the 
renormalized quark self-energy S^(fc) of Fig. [5] as 



l - z °=-{^)^ M '-Y\^^r^) = v'h ""/(# K + ^ + a-^;r (21> 

where uq is the quark spinor (uquq = 1). Because Zq is interpreted as the probability to find a bare constituent 
quark without the pion cloud, Eq. (|20[) indicates that the elementary fragmentation function is normalized to the 
number of pions per quark. This is expected from our general discussions in Section |TT] and will be elucidated further 
below. Because typical values of Zq in models based on constituent quarks are between 0.8 and 0.9, we see from 
Eq. (|20[) that the momentum sum rule dz z ^ r d q (z) will be much smaller than typical empirical values. For 
example, the NLO analysis of Ref. [ll[ found a momentum sum of ~ 0.74. From this we can anticipate that the 
elementary fragmentation functions, dg, will be very small compared to the empirical ones (see Section fV|) . 

In order to confirm that this does not mean that momentum conservation is violated, we also give the expressions 
for the distribution function of a quark q inside a parent quark Q and for the fragmentation function of q — > Q. 
The operator definitions of these functions [fS{x) and D®{z)) are exactly the same as in Eqs.(l) and (2) with the 
replacement h — > Q, where the state \p(Q)) refers to fixed flavour, spin and color (c.f. the comments in footnote [3]). 
Again we will use the symbol to denote the elementary fragmentation process. The relevant cut diagrams are 
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p-k 

= > • > X > • > 




FIG. 3: Figure (a) depicts the cut diagram (left) and Feynman diagram (right) for the loop term in f®(x) of Eq. (122f> . Here 
k- = xp- and the two quark lines with momentum k are connected by a 7 + . Figure (b) depicts the cut diagram for the loop 
term in d®(z) of Eq. (|23[) . Here p_ — zk- and the two quark lines with momentum k are connected by a 7 + . This diagram 
refers to a frame where p T — and the substitution in Eq. (0 is performed in the final transverse momentum integral. 



shown in Fig. [3] and a straightforward calculation, following the rules already indicated in Eqs. (|13[) and (fT5)) , gives 8 



mx) = Z Q 5(x-l)5 q , Q +[ - 



1 r q T Q \ 3 2/1 ^ f d 2 k T k 2 T + M 2 (l-x) 2 



(2^) 3 [fc| + M 2 (l 



g^Q<5(z-l)(5 9 ,Q + 



1/1 T q T Q ^ 3 „2 



£(1-*) 



,2 6 7 2 
In accordance with Eq. © these relations show that 



pi + M 2 (l-z) 5 



(2tt) 3 [p2 + M 2(! _ z )2 + z m 2 
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(22) 
(23) 

(24) 



Therefore the two quantities in Eqs. (|22j) and (f2"3")l describe essentially the same object, namely the splitting function 
of a quark to another quark, which also includes a "non-splitting" term proportional to Zq. The normalization is 



z) = Z Q + (l- Z Q ) = 1, 



(25) 



where the factor 6 represents the summation over the spin and color of Q. As expected, the second term in Eq. 
can be obtained from the elementary q — > it fragmentation function expressed in Eq. (fTT)) , via the substitutions 
z — > 1 — z and r w — > (r q — tq)/2. This directly leads to momentum conservation for the fragmentation of q into either 
Q or vr (see Eq. (JUJ)). 

This connection between splitting functions can also be viewed another way: The second term in Eq. (|22p . which 
describes the distribution of q inside Q with a pion spectator, suggests that via the substitutions r q /2 — ► tq/2 — 
and x — > 1 — a; we obtain the distribution function of a pion inside the quark Q, namely 



„ Q . , 1 , 2 f d 2 k T k 2 T + M 2 x 2 

2 •/ ( 27r ) fe| + A/ 2 :r 2 + 1 



x) m 2 ] 



(26) 



The tree level terms proportional to Zq in Eqs. i ]22[ ) and (123 I t come from the vacuum state in the sum over n in Eqs. JTJ and J2]l, which 
contributes for the case where p(h) is a quark. Using tp = •/ Zqj/j, where j/j is the renormalized quark field with unit pole residue of the 
propagator, gives the Zq terms in Eqs. {22) and J23j- Note, in the loop terms all factors Zq of the propagators cancel. We also note 
that the loop terms in fq{x) and (z) formally satisfy the DLY relation, that is d® loop (z) = (— z /6)f®\ oop ( x = however it is 

violated after regularization. 
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Comparison with Eq. (fT"7f gives d^(z) = f%(z), in accordance with Eq. ([8]). This relation further elucidates the 
interpretation of the normalization given in Eq. (|20[) as the number of pions per quark, namely 



[ dz Y, W = I dz = i - Za 

Jo _ Jo 



(27) 



Finally, we write down the momentum sum rules for the elementary splitting functions. In terms of the distribution 
functions we have 



[ dxx \Y,f?(x)+J2f®(x) ) =Z Q + [ dxxJ2f®(l-x)+ [ dxxJ2f®(x) = l, 

J ° \r q r „ J J ° r, ^0 



(28) 



where in the second equality we used x — > 1 — x and Eq. (|27| . In terms of the fragmentation functions Eq. (|28p 
becomes 

\zz [6 + =1 - (29) 



In reference to the form of Eq. (|23[) . we have the following simple interpretation of the momentum sum rule of Eq. 
Because Zq is the probability that the initial quark q does not fragment at all, the fraction Zq of the momentum 
stays with the initial quark. The remaining fraction (1 — Zq) is shared among the quark remainder and the produced 
pion, that is, the first and second terms in Eq. (|29|) . 

Although a description of fragmentation functions using only the elementary fragmentation processes does not 
violate any conservation law, it is completely inadequate for the following reasons: Firstly, there is a large probability 
{Zq) that the initial quark does not fragment. Secondly, if it does fragment the momentum fraction 1 — Zq is shared 
between the quark remainder and the pion. Both points are in contradiction to the usual assumption of complete 
hadronization, which is expressed by the momentum sum rule of Eq. (|10[) . 



IV. GENERALIZED PRODUCT ANSATZ FOR QUARK CASCADES 

From the previous section, it is clear that we have to consider the possibility that the fragmenting quark produces a 
cascade of mesons. A simple model to describe cascades is the quark jet-model of Field and Feynman [l9(. However, the 
product ansatz used in this model assumes from the outset that the probability for fragmentation in each elementary 
process is 100%, and that the quark produces an infinite number of mesons. Because these assumptions are inconsistent 
with our present effective quark theory, we will first introduce a generalized product ansatz, then explain its physical 
significance and its relation to the original quark jet-model. 

We assume that the maximum number of mesons which can be produced by the fragmenting quark is N . We 
then consider a process where the initial quark with light-cone momentum fc_ = Wo (which we will simply call the 
momentum in the following) goes through a sequence of momenta Wo ^ W\ ^ W2 ^ • ■ • ^ Wn, and introduce the 
momentum ratios 

^ = ^F^> n=l,...N. (30) 

Our product ansatz for the fragmentation function, which we will motivate shortly, is: 

D K Z ) = Y J / d Vi / dm--- d VN 
m=1 Jo Jo Jo 

Qdf{ m ) -Qd Q Q l{ m ) 6^(^)5^ - z rn ) 5(t«, (tq^ - r Qm ) 12) . (31) 

Qn 

Here the functions dg (77) are our elementary Q — > Q' splitting functions of Eq. ([23]) , which represent the probability 
that a quark of flavour Q makes a transition to the quark Q' , leaving the momentum fraction 77 to Q 1 . A sum over 
repeated flavour indices is implied in Eq. (|3ip : a flavour sum over the quark remainder (Qn) is included; for the case 
N = 1 we define Qo = Q', and the symbol S(i,j) denotes the Kronecker delta. The factor 6 which multiplies each 
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-M>->-o- 

W Wi 



W,r, 



WoZ 



W m 



Woz 



Wn 



k 



Wo Wi 



Wm-l W m 



w k 



> o > ZZ > X > + > * > 

Zq 

FIG. 4: The left hand side of the top figure is a graphical representation of Eq. (|3ip and the right hand side of this figure 
represents Eq. (|37[1 . The open circles denote the elementary q — » Q fragmentation function of Eq. (|33[) and the dots represent 
the second (meson emission) term in Eq. (|33p . In the m th step, where a meson with momentum zWo is selected by the 
delta-function in Eq. (|31[) , only the meson emission term contributes. The term P(k) is the binomial distribution of Eq. (|38[) 
and the squares represent the renormalized meson emission term, F®(z), given by Eq. (140 \ . The bottom figure is a graphical 
representaion of Eq. (|33l) . 



elementary splitting function comes from the sum over spin and color. The delta function in Eq. (|3Tj) selects a meson, 
which is produced in the m th step with momentum fraction z m of the initial quark: 

Wm-l — Wm n \ i 1 i i 

z m = — =??;l-772 »7m-i • (1 - Vm), where m > 1, and zi=l-r?i. (32) 

Because the pion has a mass we will exclude the unphysical case of z = 0, that is, whenever a pion is produced in the 
m th step we will assume that r\ m 7^ 1 in Eq. (f32|) . 

We will write the q — » Q splitting function of Eq. (|23l) . including the spin-color factor 6, in the form 

6^(z) = Z Q «S(z-1)5,, q + F«(z), (33) 

where 




The function F satisfies the normalization (see Eq. (|25|) ) 

[ dzF^(z) = [ dzF(z) = 1 - Z Q . (35) 

q JO JO 

For the case N = 1 it is easy to see that Eq. (|3"Tj) reduces to the elementary fragmentation function of Eq. (j 1 T[) . 
namely 

D*{z) ^ F«(l - z)| TQ=T? _ 2r „ = 1 (1 + r,r ff )F(l - z) = dj(«). (36) 

In order to illustrate the physical content of the ansatz expressed by Eq. ([3~T|) we rewrite it identically as follows: 
Noting that each factor of the product in Eq. (f3"Tj) consists of the two terms in Eq. (|3"3"j) . it is easy to see that all products 
with the same number (call it k) of F's and (N — k) number of Zqs make the same contribution to D^(z). Therefore 
we can introduce an ordering of the 77's in Eq. (|3"Tj) . That is, take the first k 77's not equal to one (771, 772, • • • Vk 7^ 1), 
and the remaining 77's equal to one (77/0+1, Vk+2, ■■■Wn = l)j multiply by the combinatoric factor f 1 ?) and perform a 
sum over k. For some fixed k, only terms with m ^ k will contribute to the sum in Eq. (pTTj) , because z m of Eq. (|3"2"]) 
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must be non-zero. 9 Then Eq. (|31[) is rewritten identically as 



jv jv i i i 

W = E E p ( fc ) / d K / ^ . . . / 



N 



dr)k 

E P ? X M F^{ m ) . . . S(z - z m )S(r n , (tq^ - r Qm ) /2) , 

= E d ;,hW. ( 3? ) 

m— 1 

which is expressed graphically in Fig [?1 The binomial distribution 

P(k)= (T\z$-"(l-Z Q ) k , (38) 
is the probability of producing k mesons out of a maximum of N mesons and satisfies the normalization condition 

N 



E p ( fc ) = L ( 39 ) 



fe=0 



In Eq. ((371) we defined the renormalized function = F®/(1 — Zq), that is (see Eqs. d34j) and p5)l ) 

F?^)=(\- TJ ?]H^ where F{z) = and (40) 



2 6 



[ dzJ2 F ?(z)= [ dzF(z) = l. 

JO r, JO 



(41) 



The physical interpretation of Eq. (|37|) is as follows: 

• P{k) is the probability that k mesons out of a maximum of N mesons are produced. 

• Fq (77) is the probability density that, if a meson is emitted from the quark Q, the momentum fraction ij is left 
to the remaining quark Q' . 

• The product F(j)i) -F (772) ■ ■ ■ F{rjk) is the probability density that, if k mesons are produced, each meson carries 
its momentum fraction z m (m = 1, . . . k) of the original quark, where z m is given by Eq. (|32p . 

• / m \(z) is the probability density that the m th meson has the momentum fraction z of the original quark. 
This implies that at least m mesons must be produced, which corresponds to the lower limit (k — m) of the 
summation in Eq. (|3T|) . The total fragmentation function D^(z) is then obtained by summing the probability 
densities D* q (m) (z). 

We note that the original ansatz of Field and Feynman 19] is an infinite product, which formally emerges from 
Eq. (|3"T|) if we take the limit N — > 00 and assume that P(k) is equal to zero for any finite k, that is, the probability 
of the fragmenting quark to emit a finite number of mesons is zero. 

We now proceed with Eq. (|37[) in order to find the integral equation satisfied by the fragmentation function. For a 
fixed m, we can integrate over ?y m +i, . . . tjn by using the normalization of F, that is, 



Jo n Jo Jo n , 



(42) 



As explained earlier, we only consider the case z > 0. 
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Then for all k ^ m the integrations over the same variables rji, . . . r\ m remain, and the sum over k refers only to the 
probabilities P(k). Performing the shift rj m — > 1 — rj m in the integral over ?y m , we obtain 



Dl (m) {z) 



VP(fe) j f dm ['dm... f 

,k= m ) Jo Ja Jo 



dr]„ 



f? 1 (m)F^ (m) ■ ■ ■ F%zl fam-i) d Qm _ x (vm) S(z - mm ■ ■ ■ r? ro ). (43) 

The function d^(z) = c?J(z)/(l — Zq) is the renormalized elementary q — > 7r fragmentation function, therefore (see 
Eq. pU) 



<£(z) = - z)\ TQ=Tq ^ = I (1 + ^)^(1 - z). 

From Eq. (|43| it is easy to derive the following recursion relation for m > 1: 



^ (n 0(*) = -Rr. 



^ Q ^Q( m -l) 



(z), where m > 1, 



while for m — 1 we have 



We have introduced the following ratios: 



where n = 1, 2, . . . N, 



and used the following notation for the convolution of two functions A(z) and B(z): 



[A®B](z) = J dzi I dz 2 5{z- z 1 z 2 )A(z 1 )B(z 2 ). 
The total fragmentation function then becomes 



D^{z) = R 1 d^(z) + Y,R r> 



n=1 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



where D^,-* can be obtained from the recursion relation of Eq. (|45l , with the starting value given by Eq. (|46|) . 

It is interesting at this stage to derive the sum rules for the fragmentation function. A simple calculation using 
Eq. (|43j) gives the following expressions for the multiplicity, the momentum sum and the isospin sum: 



-1 iV 

/ dz ]T D" q {z) = kP(k) = N(l - Zq), 

Jo r„ fc=l 

r 1 N / 

/ dzY,zD*(z) = l-J2P(k)(zF) k = l- (z Q + (l-Z Q )(zF) 
Ja r, ' k =o 

r 1 I N f ^ 



1 - ( Z Q - -(1 - Zq) 



N 



(50) 
(51) 
(52) 



where (A) = dzA[z). These expressions can be understood as follows: If k mesons are produced with probability 
P(k), then Eq. (|50|) is simply the mean number of mesons; the quantity P(k)(zF) k in Eq. (|5ip is the mean momentum 
fraction left to the quark remainder; and the quantity P{k) (—1/3) in Eq. (|52[) is the mean isospin fraction left to 
the quark remainder. 

Eqs. (|51J) and (|52p indicate that, in the present model, it is not possible to transfer the total momentum and isospin 
of the original quark to the mesons, if the maximum number of mesons is finite. The momentum and isospin sum 
rules given in Eqs. (|10p and (jlip are valid only in the limit N — > oo. While this may indicate a conceptual limitation 
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of the jet-model, we note that in general, the QCD based empirical analysis of fragmentation functions also leads to 
divergent multiplicities. Therefore, we find it more important to satisfy the momentum and isospin sum rules given in 
Eqs. (fT0|) and (fTT|) than to have finite multiplicities, and therefore we take the limit N — > oo. The results then become 
independent of the form of the distribution P(k), if the following condition is satisfied for the ratios in Eq. (|4"T|) : 

R n JV =5° 1, for all n = 1, 2, . . . (53) 



In fact, it is well known that in the limit N — > oo our binomial distribution of Eq. (|38p becomes a normalized 

_(h-k ) 2 

Gaussian distribution (normal distribution) ^ 2 = e ^ , with the same mean value ko = N(l — Zq) and variance 

c 2 = NZq(1 — Zq) as the original binomial distribution. The validity of Eq. (|53"|) can then easily be confirmed. In 
fact, any distribution which approaches a normal distribution in the limit N — > oo satisfies the condition given in 
Eq. (MD. 10 

Using Eq. (|53|) . wc sec from Eq. (l49|) that our fragmentation function satisfies essentially the same integral equation 
as in the original quark jet-model [19j : 

DZ(z)=d*(z)+[F?®D Q ](z), (54) 

where the driving term is given by Eq. (|44[) and the integral kernel by Eq. (|40[) . We finally write down the equations 
which we solve in the next section. Defining two functions A(z) and B(z) by the isospin decomposition 

D^(z) = ^[A(z)+T q r n B(z)}, (55) 
and using Eqs. (|4"0"|) and (|4"4"|) . we find the following integral equations for A(z) and B(z) from Eq. ([54")) : 



A(z)=F(l-z)+ / — F I - ) A(y), (56) 

j z y \y ; 

B(z)=F(l-z)-l t^-Fl- ) Bun. C7) 



where F(z) is obtained by renormalizing the function F(z) in Eq. ([34]) to unity. Using Eq. (|55p . we obtain the following 
expressions for the favoured, unfavoured and neutral fragmentation functions: 

D? =Dl=Dl =Df = ±(A + B), (58) 
Dl = Df = Df =Df = ±(A-B), (59) 
DZ° =Df =Df =Df = ±A. (60) 

From the form of Eqs. ([56|) and ([57j) it is easily seen that (zi) = 1 and (B) =3/4, which leads to the momentum and 
isospin sum rules of Eqs. (fT0|) and (fTTjl . For large z, both functions A(z) and B(z) approach F(l — z) and therefore 
the unfavored fragmentation functions in Eq. ([5^]) are suppressed for large pion momenta. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



In this section we present the numerical results for the fragmentation function of Eq. (I54|) in the NJL-jet model. 
For reference, we also give the results for the elementary distribution function of Eq. (|14j) . Because the application 
of the NJL model to the calculation of the quark distribution functions in the pion has been explained in detail in 



The fact that in the limit N — » oo the binomial distribution becomes a normal distribution is known as the Moivre-Laplace theorem, 
which can be formulated rigorously in integral form ( "weak convergence" ) . The central limit theorem [3l] | is an extension of the Moivre- 
Laplace theorem to general distributions P(k) with mean value proportional to N and variance c 2 oc N. This indicates that Eq. H53H is 
actually valid for a wide class of distributions. Although our NJL-jet model ansatz of Eq. I|31|) leads to the binomial distribution, in the 
limit N — » oo the results hold for a wide class of distributions. 
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Ref. [22j, we will not repeat the explanations of the model here. For convenience, we will use the same regularization 
scheme, namely the invariant mass, or Lepage-Brodsky (LB) [32| regularization scheme, with the same parameters as 
in Ref. [22j ■ The LB scheme is suitable for regularizing integrals in terms of light-cone variables and in terms of the 
usual variables it is equivalent to the familiar 3- momentum cut-off scheme [22j ]. That is, if we denote the 3-momentum 
cut-off by A3, which is fixed in the usual way by reproducing the experimental pion decay constant, a bubble- type 
loop integral with two intermediate particles of mass M\ and M2 is regularized by cutting off their invariant mass 
M12 according to 



In terms of light-cone variables, if we associate with particle 1 the transverse momentum q T and the momentum 
fraction y of the total P_ momentum, and to particle 2 we associate the momentum fraction (1 — y) and transverse 
momentum — q Tl then their invariant mass squared is 



The requirement Mi 2 ^ A12 then leads to a y-dependent transverse cut-off: q% ^ A 2 2 y(l — y) — M 2 (1 — y) — M^y. 
This condition also restricts the values of y from below and above (0 < y\ ^ y ^ 2/2 < 1)- For example, for the 
integral in Eq. (|17|) of the elementary q — > 7r fragmentation function we have Mi = and M2 = M, for the integral 
in Eq. ()23|1 of the elementary q — > Q fragmentation function we have Mi = M and M2 = and for the integral in 
Eq. (TT3]) of the distribution function we have Mi = M2 = M. We also note that this regularization scheme does not 
violate the sum rules. 

Following Ref. [22j we use a constituent quark mass of M = 300 MeV. Then A3 = 670 MeV and the invariant mass 
cut-offs for the (7r,q) and (q,q) systems are 1.42 GeV and 1.47 GeV, respectively. We did not investigate whether 
other parameter sets or other regularization schemes lead to a better description of the fragmentation functions. 

As usual, we will associate a low energy renormalization scale (Qq) to our NJL results and evolve them in Q 2 by 
using the QCD evolution equations. For the evolution of the fragmentation functions we limit ourselves to LO. In 
this case it has been verified [ItJ that a formal application of the DLY relation, see Eq. (|12jl . leads to the correct 
connection between the evolution kernels of the distribution and fragmentation functions (see Appendix B). However, 
the DLY relation is not actually used to relate the distribution and fragmentation functions themselves. We therefore 
use the Q 2 evolution code of Ref. [33| at LO for the distribution functions, and perform the transformation of the 
kernels as explained in Appendix B to obtain the LO evolution of the fragmentation functions. 11 

In Fig. [5^ we recapitulate the results of Fig. 4 of Ref. [22], and show the minus-type (valence, q — q) it-quark 
distribution in a tt + and in Fig. [SJd we give the result for the plus-type (q + q) it-quark distribution in a tt + . The 
dotted line shows the NJL model result based on Eq. (| 14(1 . the solid lines illustrate the distribution obtained by 
associating a low energy scale of Qq = 0.18 GeV 2 to the NJL result and performing the Q 2 evolution at LO and NLO 
to Q 2 — 4 GeV 2 . The dashed line shows the empirical NLO parametrizations of Ref. |8|. We see that the LO and 
NLO results show quantitative differences because of the rather low value assumed for Q 2 , although the qualitative 
behaviours are similar. 

In Figs. Owe present the corresponding results for the minus- type and plus- type fragmentation functions for u — > tt + . 
The NJL-jet result, given by the dotted line, is the solution of the integral equation in Eq. (|54| . Therefore the dotted 
line in Figs.EK and[6}} show the functions jB(z) and ^A(z), respectively (see Eqs. (|58|) and (|59|) ). In order to see the 

importance of the cascade processes, we also plot the driving term of the integral equation, namely I F(l — z), as the 
upper dash-dotted line, which is the renormalized elementary fragmentation function. As the lower dash-dotted line 
we illustrate the elementary fragmentation function, namely f -F(l — z). The result of the evolution of the dotted line 
(Qq = 0.18 GeV 2 ) to Q 2 — 4 GeV 2 at LO is shown by the solid line and the dashed line shows the empirical NLO 
result of Ref. llj, evolved to Q 2 = 4 GeV 2 . 

Several important points are illustrated in Figs. [5] Firstly, as anticipated in Section IIIH the elementary frag- 
mentation function (lower dash-dotted line) is very small. Secondly, Fig. [5Jd shows the tremendous enhancement at 
intermediate and small z of the plus-type fragmentation function caused by the cascade processes (iterations of the 
integral equation of Eq. (|54p ). while for the minus-type fragmentation function of Fig. [6^ a small reduction is seen. 




(61) 




M 2 + q 2 T M 2 + q 2 T 

y 1 - v 



(62) 



11 The DLY based relation between the evolution kernels for distribution and fragmentation functions is violated at NLO |17| . Unfortu- 
nately, a NLO evolution code for the fragmentation functions is not yet publicly available. In this paper we do not attempt a quantitative 
comparison with the empirical functions, therefore we leave the NLO calculation for future work. 
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Q 2 =4 GeV 2 




NJL (Ql = 0.18 GeV 2 ) 
NJL 
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0.4 




FIG. 5: Figure (a) depicts the minus-type (valence) quark distribution x{fZ (%) — fS {x)) and figure (b) illustrates the 
plus-type quark distribution x{fZ {x) + f£ (x)) of the u-quark in a ir + . The dotted line is the NJL model result, used as input 
{Ql = 0.18 GeV 2 ) for the Q 2 evolution. The solid line labeled by LO (NLO) is the result of LO (NLO) evolution to Q 2 = 4 
GeV 2 . The dashed line is the empirical NLO result of Ref. @ at Q 2 = 4 GeV 2 . 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 



(a) z (b) 

FIG. 6: Figure (a) depicts the minus-type fragmentation function z{DZ + (z) — -DJ + (z)) and figure (b) illustrates the plus-type 
fragmentation function z{D™ + (z) + (z)) for u — > n + . The dotted line is the NJL-jet model result, used as input (Q 2 , = 0.18 
GeV 2 ) for the Q 2 evolution. The lower dash-dotted line is the elementary fragmentation function (d, of Eq. (|17p 1 and the 
upper dash-dotted line is the renormalized elementary fragmentation function (d£ of Eq. (|44p ). which is the driving term of 
the integral equation expressed in Eq. (|54p . The solid line is the result after LO evolution to Q 2 = 4 GeV 2 . The dashed line is 
the empirical NLO result of Ref. [TO], evolved to Q 2 = 4 GeV 2 . 



Thirdly, the calculated result shown by the solid line has the correct order of magnitude for intermediate and large z, 
when compared with the empirical function. This point, which reflects the fact that our model satisfies the momen- 
tum sum rule, is very important, because effective quark model calculations completed hitherto only considered the 
elementary fragmentation functions and introduced some ad hoc parameters (like normalization constants) to obtain 
the correct order of magnitude. Quantitatively, Figs. [6] indicate that our fragmentation functions are too big at large 
z and too small at smaller z. This is natural for the following reasons: Firstly, we can expect that a NLO calculation 
will lead to a softening of the fragmentation functions. Secondly, some of the observed pions are secondary ones, 
which come from the decay of primary p and u> mesons. Thirdly, the coupling to other fragmentation channels, in 
particular the nucleon, antinucleon and kaon, will transfer some amount of the hard quark momentum to these other 
hadrons. Also, one should not forget that the empirical fragmentation functions have very large uncertainties, which 
are not indicated in our figures. Nevertheless, Figs. [5] indicate that the present NJL-jet model provides a reasonable 
starting point for the description of fragmentation functions. 

Fig. [7k. shows the results for the favoured fragmentation function of Eq. ([55)) and Fig. [7b shows the unfavoured 
fragmentation function of Eq. ([55)) . Note, these figures correspond to half the sum and half the difference of the 
curves in Figs. [5] The upper dash-dotted line in Fig. [7k shows the driving term, §-F(l — z), of the integral equation in 
Eq. ([54]) . and the lower dash-dotted line shows the elementary fragmentation function, §-F(l — z). For the unfavoured 
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FIG. 7: Figure (a) depicts the favoured fragmentation function zDZ (z) and the figure (b) illustrates the unfavoured frag- 
mentation function zD£ (z). In figure (a) the lower dash-dotted line is the elementary fragmentation function (dj of Eq. (|17[) ) 
and the upper dash-dotted line is the renormalized elementary fragmentation function (dj of Eq. <|44[) ). which is the driving 
term of the integral equation in Eq. (|54[) . Note, these functions are zero for the unfavoured case. The solid line is the result 
after LO evolution to Q 2 = 4 GeV 2 . The dashed line is the empirical NLO result of Ref. [TH, evolved to Q 2 = 4 GeV 2 . 



case these two functions are zero. Both figures demonstrate the importance of cascade processes in the present NJL-jet 
model. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we used the NJL model as an effective quark theory to study the simplest fragmentation function, 
namely, the fragmentation of unpolarized quarks to pions. Our aim was to develop a framework which satisfies the 
momentum and isospin sum rules in a natural way, without the introduction of ad hoc parameters. This framework 
should also give fragmentation functions that have the correct order of magnitude at intermediate and large z. We 
explained in detail, that for this purpose, the simplest approximation where a truncation is made to the one-quark 
spectator state, in the denning relation given by Eq. (O, is completely inadequate. Although this approximation 
does not violate any conservation law, it gives very small fragmentation functions; because the probability for the 
elementary fragmentation process is small in effective theories based on constituent quarks and the quark remainder 
can carry an appreciable amount of momentum. 

In order to overcome these difficulties we followed the idea of the quark jet-model and made a generalized product 
ansatz to describe the cascade processes in the NJL model. We explained that this ansatz corresponds to a binomial 
distribution for the number of mesons emitted from the quark. However, in the limit that the maximum number of 
mesons becomes very large the results are independent of the form of this distribution function. Our formulation thus 
represents an extension of the original quark jet-model, which assumed an infinite number of mesons from the outset. 
We have shown in detail that this NJL-jet model describes fragmentation processes where 100% of the initial quark 
light-cone momentum is transferred to mesons. The momentum sum rule of Eq. (|10p . which is assumed valid in all 
QCD based empirical fits, is then satisfied automatically without introducing any new parameters into the theory. 
We have also shown that the isospin sum rule of Eq. (|TT|) is naturally satisfied in this approach. 

The comparison with the empirical fragmentation functions shows that our calculated functions have the correct 
order of magnitude for intermediate and large z. We highlighted that a straightforward extension to include the 
NLO terms in the Q 2 evolution and to include the effect of primary p and uj mesons, as well as fragmentation to 
other hadronic channels, will improve the description. Therefore we can conclude that our NJL-jet model provides a 
reasonable framework to analyse fragmentation functions in an effective quark theory. 

For future work in this direction it is important to derive the jet model type product ansatz from field theory. 
The rainbow-ladder scheme for the quark self-energy may provide a suitable framework for this purpose. An attempt 
can then be made to use this truncation scheme to consistently describe the cascade processes for the fragmentation 
functions and to include the contribution from the hadron cloud around the quark for the distribution functions. 
However, it is important to bear in mind that a truncation scheme which works well for fragmentation processes may 
not be suitable for the distribution functions and vice versa. To establish a scheme which respects the sum rules and 
which gives a satisfactory description of both types of processes is an important task for future research. 
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APPENDIX A: PROOF OF THE DLY RELATION 

In this Appendix we will prove the DLY relation expressed in Eq. (Ti"2"|) in two independent ways. First, we follow 
the original derivation of Ref. [l6l ] in terms of the hadronic tensors and second we start from the operator definitions 
given in Eqs. (HJ) and @. In order to illustrate the spinor algebra the formulae in this Appendix refer to the case 
where h is a proton, however it is trivial to modify the expressions for the case where h is a pion. 

1. General crossing relations 

We consider the following Green function 

Ml(p,Pn) = J d 4 xe-^{ Pn \T(O a (O)^ (x)) |0), (Al) 

where $p(x) is an interpolating field for the nucleon and O a is another local field operator. We also define the 
iV-amputated Green function by Mp(p,p n ) = T (p,p n ) iGN^pip), where Gn is the nucleon propagator. From the 
spectral representation of Eq. (|A1|) or from the familiar reduction formalism, we can derive the relations 

( Pn \O a \p) =T a (p, Pn )y/2M^u N {p S ), (A2) 

(p,p n \O a \0) = {±)T{-p,p n )^/2M^v N (ps). (A3) 

In Eq. (|A3[) the sign is (+) if O is a fermion type operator and (— ) if it is a boson type operator. Also, p denotes 
an antinucleon with 4-momentum p^ = (_Ejy (p),p). The nucleon spinors are denoted by ujv and v/v- Our covariant 
normalization implies the following matrix elements of the nucleon field operator: (0|$(0)|p) = v / 2Mjv un{ps) and 
(p|$(0)|0) = \/2Mjv vn(ps). Eqs. (|A2[) and (|A3|) are the basic crossing relations which will be used in the following. 

2. Comparison of hadronic tensors 

Here we use the above crossing relations to find the connection between the hadronic tensors (spin-independent 
parts only) for the processes eh — > e'X and e + e~ — > hX, where h denotes a hadron (proton) Q: 



wT(p>?) = ^£ n (2n) 4 6 4 (q+p-p n )(p\j»\p n )(p n \r\p) 

= + q ^f) 9 2 ) + jr q {f- ^V) {f - ^V) F ^ x > (A 4 ) 

WT(p J? ) = (^) 4 S\q-p-p n )(0\J>,P^)(P^\J V \0) 

= (-»"" + ^ *<•■«"> + JTq if - IF*") - 7^) ,A5) 
Here \p) is the state of the hadron h with momentum p and we use x — 7^- and z = — __I We also defined 

J2 n = J2 n I (2n) s ^^Pn ~ ^n) ®{Pno)i where M n is the invariant mass of the state n. Using Eq. (|A2|) and its complex 
conjugate for the current operator J v : 



{p n \r\p) = V2M^T (p ,p n )u N (ps), (A6) 



(p|J>„) = ^2M^u N {ps)Y v {p,p n ), (A7) 
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where — (70^")^, that is, — V^j . We insert these relations into Eq. (|A4p . Since we consider the spin- 
independent part only, we can sum over the nucleon spin s and divide by 2, using ^ s un(ps)un(ps) — ^m^T • This 
gives 

47rWT(P.?) = ^E {^f5% + P - Pn )TT \(p + M N )T^p,p n )r(p, Pn )] . (A8) 



For the hadronic tensor in Eq. (|A5[) . we first use charge conjugation and then Eq. (|A3|) and its complex conjugate for 
the current operator J M : 



(0\J»\p,K) = (OIC- 1 (CJ^C-^Clp,^) = (0| (CJ^C- 1 ) \p, Pn ),= -(0\J»\p, Pn ) = y/2M^v N (p S )T^{~p,p n ), (A9) 
(p^J^O) = -<p,p n |J M |0) = v / 2Mw ; r /i (-p,}>„)i;^(iM). (A10) 
We insert these relations into Eq. (|A5[) . sum over the nucleon spin s and divide by 2 using ^ s vn(ps)vn(ps) — 

- ZJ ^mr- This s ives 

Wr(p,9) = -\j2 n (2ir) 4 S 4 (q-p- Pn )Tr [{-f + M N )T»(-p,p n )T" (-p, Pn )\ . (All) 

By comparing Eqs. (jA8|) with (jAllll we obtain the DLY crossing relation for the hadronic tensors: 

W^(p,q) = -W^(- P ,q), where s h = ±. (A12) 

The minus sign in Eq. (|A12|) comes from the Dirac algebra, and for a spinless hadron the minus sign is changed to 
plus. Eq. (|Al2j) implies the following relation between the structure functions in Eqs. (|A4j) and (|A5j) : 12 

F'Kz, q 2 ) = q 2 ) = -If ( V) > (A13) 

T' 2 \z, q 2 ) = F>\-x, q 2 ) = F 2 h (I, q 2 ^j . (A14) 
The well known relation F 2 (x) — 2xF[ l {x) becomes, with x — > — x and using the first equalities in Eqs. (|A13|) and 

F 2 h (z) = --F^(z), (A15) 

z 

which also holds for spinless bosons. 

The connection between the structure function F 1 and the fragmentation function D^(z) in the Bjorken limit is as 
follows: The cross section for the process e + e~ — > /iX is Q| 13 

do-' 1 2a 2 7T2; /— fc. 2 — h. ,,\ 4 a 2 irz — h, , , , „ „. 

"d7 = -£3- [ F i (*> 9 ) + e ^2 9 2 )) = 3 — 5- *i (*)• (A16) 

Usually this is divided by the total cross section for e + e~ — > hadrons 

47ra 2 v-^ , 4-7ra 2 



Ctot 



E^^ (AIT) 



2 ^-8- 3 
9 



By relations like Eq. (I A13t we mean the following: Take a particular physical value of 2 for the (e+, e ) process (0 < z < 1). Then the 
corresponding (unphysical) value of the Bjorken variable for the (e, e') process is x = 1/z and Eq. <Xl3t gives the connection between 
the structure functions. 

We remind the reader that the symbol h denotes a particular hadron with a specified spin direction, e.g., p f (although the spin averaged 
cross section considered here does not depend on the spin direction). Therefore, the cross section measured for the case that the spin 
of the produced nucleon is not observed has an additional factor of 2, which is not included in Eq. (I A16E . 
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where ^ refers to the quark flavour only. Then we obtain 



otot dz R 

This is compared to the original definition of the fragmentation function [l[ 

1 da h 
ctot dz R 

to obtain 



1 da h 1 —h, . 

— zF 1 (z). (A18) 



UY,4 {D h q {z)+ Dh M ( A1Q ) 



F i( z ) = - z J24 (° h M + Dh M • ( A2 °) 
i 

Because we know how to express F^ix) in the Bjorken limit by the distribution functions f^{x), we obtain from 
Eq. JXnH :: 

^w-f (i)-5?«!(jf(i) + Jf(i))- < A21) 

Comparing (|A20[) and (|A21|) we obtain 

D fc) = -ftf(:;)> where Sft = 5' (A22) 

and a similar result holds for the antiquarks. Eq. (|A22[) expresses the DLY relation of Eq. ([12]) between the distribution 
and fragmentation functions. For the case of a spinless hadron the minus sign in Eq. (|A22[) becomes a plus sign. 



3. Comparing the operator definitions 

Starting from the operator definitions given in Eqs. |T]) and (O, we obtain 

f«( x ) = \yi (5(j ) -a:-P-+Pn-)(p|?|Pn)7 + (Pn|^b} 1 (A23) 

For definiteness wc consider again the case where h is a proton. We use O a — ip a in Eq. (|A1[) . which gives 

ujv(ps), (A25) 
v N (ps). (A26) 

We insert Eq. (|A25j) and its complex conjugate into the operator definition, Eq. (|A23|) . and average over the nucleon 
spin. This gives 

^ = \J2J(P- X -P- +P«-) Tr W + M J v)r(p,p, i ) 7 + r(p,p n )] . (A27) 

For the fragmentation function in Eq. (|A24[) . we use the charge conjugation relations of the quark field operators 
CipaC^ 1 = (C7°) Qj g^g and Cip^^ 1 = ippCp a , where C = ij 2 j°, to rewrite the matrix elements in Eq. (|A24[) as 
follows: 

(0|Va|p,K) = (C 7 o)a/3(Pn,p|#|0)*, (A28) 
(P,KT|?a|0> = (P,Pn|#|0)C^. (A29) 

Then we use Cy C = (7^) T and Eq. ([A26]) to write 

(p,^|^ Q |0) 7 + /3 (0|^b,^) = v N (ps) [r{-p, Pn )j + T(-p, Pn )] v N {ps) ■ 2M N . (A30) 
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Averaging over the nucleon spins we finally obtain 

z 1 



D « {z) = ~h^J (l~ ~ p - ~ Pn -) Tr + M N)r(- P ,p n h + r(-p lPn )] . (A3i) 

Comparison of Eqs. (|A27|) and (|A31|) gives 



(A32) 

p^-p 



where p — > — p means to reverse all 4 components of p M and after this replacement p° = Em{p) > 0. 

We now consider the property of the distribution function in Eq. (|A27[) under p^ 1 — > — p^ . Expressing the summation 
in terms of light-cone momenta, the distribution in Eq. (|A27[) can be written in the form 



h W 4 2^J ( 27r )3 2p_(l-s) 



c>(fe + - eAr(p) + e„(p - fe)) 

5(k- - p-x) Tr \N> + M N )T(p, p - k)j + T(p, p-k)], (A33) 



where e n (p n ) = P " 2 p +M " and ejy(p) = P ^ 2 p Mjv . We then replace — * — and then fc' 1 — > — fc M in the integral. This 
gives 

<5(fc_ - p_a;) Tr + Mjv)r(-j), -p + fe) 7 +r(-p, -p + fe)] . (A34) 

Because the result of taking the trace in Eq. (|A34j) must be the plus component of a Lorentz four vector constructed 
from p^ and we have 

Tr [(-} + M N )T(-p, -p + fc) 7 +r(-p, -p + fe)] = -Tr [(|> + Mjf)T(p, p - k°h + T(p, p-k)]. (A35) 

If we use Eq. (533} to define a function by f%(x) = 6(1 - we obtain from Eqs. (fAl4|) and (535}: 

/^(iz;)|p^_p = 8(x — l)F(x). From Eq. (j A32[) we then obtain the connection between the distribution and the 
fragmentation function as 

/.» = 6(1 - x)F(x), (AM) 



(A37) 



Note, for spinless bosons there is no minus sign in Eq. (|A37j) . This result agrees with Eq. (|A22j) and would suggest 
that f£ and are essentially one and the same function, defined in different regions of the variable. 



APPENDIX B: DLY TRANSFORMATION OF EVOLUTION KERNELS 



In this Appendix we explain the DLY based relation between the evolution kernels for distribution and fragmentation 
functions, which is known to be valid at LO [ljj. Using Eq. (fl"2|) . we consider the following transformation of the 
quark and gluon distribution functions: 




where the upper (lower) sign holds if h is a boson (fermion). Using the well known evolution equations at LO [5[, it 
is easy to derive the corresponding transformation of the evolution kernels. For the minus-type (flavour non-singlet) 
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combination q — q, the kernel (Pqq) is unchanged. For the flavour singlet combination, X)i=d ft + ft; which couples to 
a gluon, the evolution kernel is transformed as follows: 



Pqq{x) Pqg(x) 
Pgq( x ) Pggi x ) 



Pqq(z) 2N f P gq (z) 



(B3) 



Here Nf = 3 is the number of flavours used in the Q 2 evolution equations. For reference, we summarize the forms of 
the individual kernels below: 



- — + -S(x-1) 

{l-x)+ 2 



4 

3 _ 

Nf [x 2 + (1 - x) 2 ] , 

4 1 + (l-a;) 2 

3 x 

x 1 — X 



6 



+ 35(1 - X) 



11 

T 



Nf_ 
3 



(5(1 -x). 



(B4) 
(B5) 
(B6) 

(B7) 
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